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I. INTRODUCTION 



During the past two decades, optimization techniques have been extensively applied to design external control 
fields to manipulate the evolution of quantum mechanical systems [U[2j[3l[4j[Sl[6l|TllSj[9l [10l [Til LL2J [13]. A 
typical application is to force the states to approach the priori prescribed targets as closely as possible, including 
both bounded [TJ ® |3J 0] and unbounded [5] situations. Also, optimization theory can be applied to improve the 
efficiency of desired quantum state transitions, e.g., the evolution time [HI [3 El El HOI 111) and the energy consumed 
by the controls \12\ 113] are most interesting. For some low dimensional quantum systems that evolve on compact 
Lie groups, one can find analytical solutions for such optimization problems with bounded pj [T^] or unbounded [5J 
controls. However, in the higher dimensional situations, numerical algorithms have to be applied. 

In this paper, we explore the optimal steering problem for the class of quantum systems whose evolution operators 
i-H can be described by the SU (1,1) matrices. The underlying system is modelled by an evolution equation of the 
£> form [2]: 

jX(t) = [A + u(t)B]X(t), X(0) = I 2 , (1) 

where X(t) is a two dimensional special pseudo-unitary matrix; u(t) is real function of time t, which is the control 
input of the system; A and B are arbitrary matrices that can be expressed as the linear combination of K X1 K y and 
t K z , which are the generators of the Lie algebra su(l, 1) and can be identified as follows: 

rS 

In |15j . Jurdjevic has taken the initial steps in the problems of optimal control for the special case when A = K x 
and B — K z (the corresponding properties of the optimal controls also can be found for the systems evolving on the 
homomorphic groups SO(2, 1) and SL(2,M.) in [T5] and [H]). Unlike the case of SU(2) [T7], however, as will be seen 
in Section [TI] system ([lj usually can't be transformed into the special form such that A = K x and B = K z . Thus, it 
is natural for us to consider the general case in detail. 

With specific realizations and representations of the Lie algebra su(l, 1) introduced, system (fTJcan be used to 
describe various of quantum dynamical processes, e.g., the superfluid system under Bose realization |l8j. the harmonic 
oscillator under ccp-realization |19j . the SU(1, 1) coherent states under irreducible unitary representation with respect 
to positive discrete series [2D]. 

The optimal control problem to be considered in this paper is formulated as follows. Given an arbitrary target 
evolution matrix Xf in SU(1, 1), we wish to find a control function u(t) that can steer the evolution matrix associated 
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with system ([I]) from its initial state I2 to some desired final state Xf, and meanwhile, minimize the quadratic cost 
function 

J(u) = [ u(t) 2 dt, (3) 
Jo 

where T is the final time. The quadratic cost index given in |3| measures the energy consumed during the steering 
process between the initial I2 and the terminal Xf. Since the remarkable difference between the quantum systems 
and the classical systems is that the evolutions of the former may be disturbed by the decoherence phenomenon. This 
is practical because that the increasing of the intensity of the electromagnetic fields, which are used to control the 
evolution of the coupled quantum system, tend to induce relaxation and decoherence phenomena. 

Based on the maximum principle for systems evolving on Lie group |21J . explicit forms of the control functions 
with respect to both normal and abnormal extremals will be derived analytically. The problem considered here can 
be viewed as the noncompact prolongation of the SU{2) case presented in [T2]. However, according to quantum 
theory and group representation theory |22j . the noncompact SU(1, 1) Lie group has only infinite dimensional unitary 
representations, and hance the associated evolution operator (or propagator) corresponding to X(t) is always infinite 
dimensional, which describes the transition between two quantum states defined in an infinite dimensional Hilbert 
space. The dynamics of the SU (1,1) coherent states is a typical example [5D]. This makes the derivation quite 
different from and far more complicated than that in the case of SU(2). 

The balance of this paper is organized as follows. In Section [TTJ some useful results including controllability 
properties and the maximum principle with respect to the quantum system evolving on Lie group SU(1, 1) are 
introduced. In Section |LU] we discuss the optimal steering problem with respect to the abnormal extremals. Properties 



of the abnormal optimal control function are characterized. In Section IV the control functions corresponding to 
the normal extremals are derived analytically for all possible cases. In Section |VJ two examples are provided for 
illustration. Finally, conclusions are drawn in Section |VI| 



II. PRELIMINARIES ON THE QUANTUM CONTROL SYSTEM ON THE LIE GROUP SU(1, 1) 

With the three generators K x , K y and K z of the Lie algebra su(l, 1) in any given evolution matrix X associated 
with system can be written as 

X = e aK 'e? K »e rK ', (4) 
where — 2n < a, 7 < 2tt, < f3 < 00. The commutation relations between K x , K y and K z are 

[K x ,K y ]=-K z , [K y ,K z ] = K x , [K z ,K x ] = K y . (5) 
With the inner product (•, •) defined by 

(M,N) = 2Tr(M7V t ), (6) 

where is the Hermitian conjugation of N, it can be verified that K x , K y and K z form an orthonormal basis of the 
Lie algebra su(l, 1). Accordingly, the drift term A and the control term B in system ([lj can be expressed by linear 
combinations of the three generators K x , K y and K z as follows: 

A = (A, K x ) K x + (A, K y ) K y + (A, K z ) K z , (7) 
B = (B,K x )K x + (B,K y )K y + (B,K z )K z . (8) 

In this paper, we say M is pseudo-orthogonal to N when (M,N^\ = 0. And M will be called elliptic (hyperbolic, 
parabolic) if (M,M^) is negative (positive, zero). Accordingly, the matrices in the Lie algebra su(l, 1) are separated 
into three different types. Since for any matrix P(zSU(l, 1), we have 

{PMP-\ (PMP-y) = 2Tr(PMP- 1 PMP- 1 ) = 2Tr (MM) = (M, Af f ) , (9) 

i.e., none of the changes of coordinates by the SU(1, 1) transformations will alter the type of an su(l, 1) matrix. Thus, 
as mentioned in Section|TJ system usually can't be transformed into the special case such that A = K x (hyperbolic) 
and B — K z (elliptic). 

We assume system ([lj is controllable in this paper, this ensures that the optimal control problem stated in the 
introduction section is always solvable. System (fil is said to be controllable on SU(1,1) if for any given target 
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evolution matrix Xf&SU(l, 1) there always exists at least one control u(t) such that X(u; T) = Xf for some time T. 
The problems of the controllability for systems evolving on Lie groups have received a great deal of attention in the 
past decades (see e.g., [23J[21]). As for the systems evolving on the noncompact Lie group SU(l, 1), a sufficient and 
necessary condition is provided in |25) and can be summarized as follows. 

Theorem II. 1 (Jgtf) 

1. If A and B are linearly dependent, then system |7J) is uncontrollable on SU(1, 1); 

2. If A and B are linearly independent, then system is controllable if and only if the set {u G R| (A + uB, A^ + 
uB^) < 0} is nonempty. 

To obtain the optimal control functions, we will make use of the well-known maximum principle on Lie 
groups [2T| [24], which, under the assumption that system (|T|) is controllable, can be summarized for the SU (1,1) 
case as follows [2"T] . 



Theorem II. 2 Assume that system |7J) is controllable, if u°(t) is an optimal control that minimizes the quadratic 
index given in |$|) and X°(t) is the corresponding optimal trajectory of |Ip. Then, there exists a constant matrix 
SGsu(l, 1) and a nonnegative real number X , not both zero, such that for almost everywhere (a.e.) the Hamiltonian 
function 

H(S; A D ; u; X a (t)) = (S, X°(t)- X [A + u(t)B]X°(t)) + l -\ u(t) 2 (10) 

is minimized with respect to u by u°(t). 

The above theory immediately provides a necessary condition for the optimality. The problems with respect to 
A 7^ are called normal, otherwise are called abnormal. 

III. OPTIMAL CONTROL FUNCTION FOR THE ABNORMAL EXTREMAL 

In this section, we will investigate the abnormal case. For this purpose, the following basic properties on su(l, 1) 
are useful. 

Lemma III.l For arbitrary pair of matrices M and N in su(l, 1), the following relations hold: 

[[M, N],M] = (M, N^) M - (M, M f ) N, (11) 
[[M, N],N] = (N, N^) M - (M, N^) N. (12) 



proof Since (11) is equivalent to (|12j because of symmetry, we only need to proof (111. If M and N are linearly 



dependent, then there exists a real constant k ^ such that M = kN. Therefore, 

[[M,N],M] = [[kN,N],kN}=0, (13) 

and 

(M, 7V f ) M ~ (M, M f ) N = (kN, N^) kN - (kN, kN^) N = 0, (14) 



i.e., the relation in (111 holds. Thus, we only need to consider the case when M and N are linearly independent, i.e., 
[M, N] ^ 0. Since 

( [M, N] i ,M) = ( [M, N] , M f ) = ( MN, M f ) — (NM, M*) = 0, (15) 

and 

(\M,N]\N) = (\M,N],N^) = (MN,N^) - (NM, 2V*) = 0, (16) 
M, N and [M, N]^ are still linearly independent. Thus we can express [[M, N], M] as 

[[M,N],M] = l i 1 M + n 2 N + ^[M,N]\ (17) 



for some constants /ii, /i2 and 113. Taking the inner product with [Af, A]t in (17 1, we obtain: 

^([M,N]\[M 7 N}^ = 0. 



Eq. (18) implys that /Lt 3 = 0, thus we can rewrite ( |17[ ) as: 

[[M,N],M] =fi x M + l x 2 N. 



Taking the inner product of (19 1 with Aft, we get 

Ml (M,Af t )+ M2 (AT, Aft) = 0. 

Since 

([[M,N],M],N^) = 2Tr([[M, AT], A/]AT) = 2Tr([M, N][M, N]) = ([M,N], [M,N?) , 



take the inner product of (19 1 with Aft, we have 

([A/, N), [Af, JV]t) = Ml (Af, ATt) + M2 (iv, iVt) 



(18) 
(19) 
(20) 
(21) 
(22) 



Make use of the equality ( [A/, N] , [Af, N] t ) = ( Af, iV+ ) - ( M, Aft ) ( AT, ATt ) (see the proof in [25] ) , from ( 20 ) and ( 22 1 
we have 

m = (M,N^), / u 2 = -(M,Mt>. 



Lemma III. 2 // Af and A^ are linearly independent and the set {u £ R\ (Af + uN, Aft + uN^ < 0} is nonempty, 
then M, N and [Af, N] form a basis in su(l, 1). 

proof Assume that [Af, N) can be linearly expressed by Af and N, i.e. there exist two real number Ai and A2 satisfy 



[Af, N] = AiM + X 2 N. 

On the one hand, making communications with Af and N respectively, one can obtain 

[[Af, N],M] = -A 2 [Af, N] = -A X A 2 M - \ 2 2 N, 
[[Af , N] , N] = Ai [Af , N] = Ai 2 M + AiA 2 A^. 



On the other hand, notice that (jl 1 1) and (12 1 would imply 



[[Af, N],M] = (Af, ATt) Af - (Af, Af f ) N, 
[[Af, N], N] = (AT, ATt) A/ - (Af, A^) N. 



Compare the coefficients of ( 24 ) and ( 25 I and those of ( 26 1 and ( 27 ) respectively, we obtain 



A x 2 = (A^Aft), A 1 A 2 --(Af,A^ t ), A 2 2 = (Af,Af t ) 



A straightforward computation, using (28), shows that, for every u € R, 

(M + uN,AP +uN 1f ) 
= (M,M^)+2u(M,N^) + u 2 (N,N^) 
= A 2 2 — 2AiA 2 u + Ai 2 u 2 
= (A 2 - Am) 2 > 0. 



(23) 



(24) 
(25) 



(26) 
(27) 

(28) 



(29) 



This contradicts with the condition that the set {u S R| (Af + uN,M^ + uN^ < 0} is nonempty. So Af, A" and 
[Af , N] are linearly independent. 

With the properties on the Lie algebra su(l, 1) obtained above, we can now draw the conclusion for the abnormal 
case as follows. 
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Theorem III. 3 Assume that the system is controllable. If (B,B^ ^ 0. then the control function u(t) — 
— (A, B^ J (B, B^ , a.e., is the only abnormal extremal for the optimization problem under consideration. Otherwise, 
there is no abnormal extremal. 



proof If X = 0, the Hamiltonian function may be rewritten as 

H(S; A ; u; X {t)) = (S, X^t)' 1 ^ + u(t)B]X°(t)) . 



(30) 



Since the minimization condition indicates that H (S; A Q ; u; X a (t)) is a.e. minimized with respect to u(t) by u (t). To 
obtain the optimal control function u (i), we differentiate H(S; A ; u; X (i)) with respect ro u and set the result equal 
to zero, which yields 



(5,l (t)^M (t)) = 0, 



By continuity, (31 1 implies that 



{S,X°(t)- l BX°(t)) = 0. 
Differentiating (32 1 with respect ro t, using ([TJ, we can obtain 

{S,X°(t)- 1 [A,B}X°(t)) =0. 



Differentiating (33 1 with respect ro t again, we have 

(S,X°(t)- 1 [[A,B] 1 A]X°(t)) + u(t)(S,X°(t)- 1 [[A 1 B],B]X°(t)) = 0. 



Utilizing (11), ( 12 1 and (31 1, we can recast (|34j) to 

{(A, Bt) + u(t) (B, Bt» (S, X^tr'AX^t)) = 0, 



a.e. 



(31) 
(32) 
(33) 
(34) 
(35) 



Notice that, according to Lemmf III. 2 A, B and [A, B] are linearly independent. If the equality (S, X°(t)~ 1 A X°(t)) = 
holds, then combining (31 ) and (33 1 we can draw conclusion that S must be zero, which is a contradiction. Therefore 
Eq.(35l implies that 



(A,B^)+u{t)(B,B^) = 0, 



a.e. 



(36) 



If (B, B^ = 0, from Eq.(36 1 we have (A, B^ = 0, which contradicts the assumption that system (fTbis controllable (see 
[2"5]). Thus, there is no abnormal extremal when (B,B^ — 0. If ^_B,_B^) ^ 0, then from ([36]) we have u(t) — 
-(A, Bt)/(B, Bt), a.e. 

Remark: The target evolution matrices that can be achieved directly by the abnormal control u = 
— (A, B^/^B, B^ are in the one dimensional Lie group corresponding to the Lie subalgebra of su(l, 1) gener- 
ated by A — (A, B^/i^B, B^B, which can never fill up the whole Lie group of SU(l, 1). Consequently, abnormal 
extremal exists only when the target evolution matrix Xf is of the form exp [c (A — (A, B^/(B, B^B)] for some 
real constant c. In order to solve the optimal steering problem subject to the terminal condition X(T) = Xf, where 
Xf is an arbitrary matrix taken from SU(1, 1), the candidates can only be normal extremals. 



IV. OPTIMAL CONTROL FUNCTION FOR THE NORMAL EXTREMAL 



In this section we explore the normal extremal control functions. According to Theorem |II.1| here we assume that 
A and B are linearly independent and meanwhile the set {u £ M| (A + uB, A^ -\-uB^ < 0} is nonempty throughout 
this section, to guarantee the controllability of the system. 

One can obtain the normal extremal, according to Theorem |II.2[ by minimizing the Hamiltonian function 
H(S; A ; u; X a (t)) as a quadratic function of u. After normalizing A D = 1, by continuity, the necessary condition 
for candidate optimal controls can be expressed as 

u(t) = — (S, X~ 1 (t)BX (t)) , a.e. (37) 

where the matrix 5* is an element in the Lie algebra su(l, 1). We introduce the following two auxiliary variables in 
the succeeding discussion: 



u A {t) = -(S,X {t)- l AX {t)), 
u c (t) ^-(S^oit)- 1 ^, B]X (t)) 



(38) 
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In order to follow standard notations, we rewrite the normal extremal u(t) in (37 1 as u B (t). Differentiate ua, u b and 
uc with respect to the time t, respectively, and make use of ([lj, ( [Tl) and (12 1, we can obtain (a.e.) 



ua = u B uc, 
u B = -uc, 



where 



u c = au A — (3u B + lu A u B — au 2 B , 



a=(A,B*), 0=(A,Ai), 7-<B,B t )- 



(39) 
(40) 
(41) 



(42) 



From (39)- (41 1, it is easy to verify the following conclusion for the normal extremal ub and the two auxiliary 



variables ua and uc- 

Theorem IV. 1 The following two quantities are conserved along the normal extremal trajectories, i.e., 

1 2 

U A + 7} U B = c i, 
1 2 1 2 

-7U4 - au A u B - - P U A = c 2 , 

for some constants C\ and c-i. 

According to Theorem IV. 1| the initial and the final values of ua, ub and uq should satisfy (a.e.) 

u A (0) + \u B {Q) 2 = u A {T) + \u B {T) 2 , 
i 7 u A (0) 2 - au A (0)u B (0) - \u c {Q) 2 -J3u A (0) 



(43) 
(44) 



where 



and 



= \lu A {T) 2 - au A (T)u B (T) - \u c {T) 2 - l3u A {T), 
u A (0) = -(S,A), u b (Q) = -(S,B), u c (0) = -(S,[A,B}), 



u A (T) = -(S,X(T)- 1 AX(T)), 
u B (T) = -(5,X(T)- 1 J BX(T)), 
uc{T) = -(S,X{T)- 1 [A 1 B]X{T)) 1 



(45) 



(46) 



(47) 



The matrix S in Eqs.(46l and (47 1 then can be viewed as parameter matrix, which has to be chosen to match the 
final condition X(T) = X f . 



Table I. The controllability of system ([lj with respect to the values of a, (3 and 



7- 



The values of a, f3, and 7 



a = 


7 < 


13 = 


uncontrollable 


(3^0 


controllable 


7 = 


uncontrollable 


7 > 


(3<0 


controllable 


j3>0 


uncontrollable 


a ^ 


7 < 


controllable 


7 > 


P < 


controllable 


/3 > 


a 2 - £7 < 


uncontrollable 


a 2 - fij > 


controllable 



the controllability of system ( 1 1 



Making use of (43) and ( |44[), from (39l-(41l, we can obtain the following differential equation for the candidate 
optimal control u(t) given in p7[) 



(u) 2 = ^u 4 + au 3 + {(3- ^c x )u 2 - 2aciu + ^c 2 - 2(3c Y - 2c 2 , a.e. 



(48) 
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where it(0) = ub(0). We will show, in the following, that the candidate optimal control function can be analytically 
solved from (48 1 in terms of the Weierstrass function. 

Since the involved system is assumed to be controllable to ensure that the optimal steering problem has solutions, 
one only need to consider the controllable situations accordingly. Table I shows the controllability properties of 
system ([lj in different cases, (see [25] for details). There are three different cases, which need to be taken into 
account, depending on the values of a and 7. 
1. Case a — and 7^0. 

In this case, the drift term A of system (JTJ is pseudo-orthogonal to the control term B while the latter is not 
parabolic. Accordingly, Eq.(48l can be simplified as (a.e.) 



(uf 



+ {0- 1Cl )u 2 + 7 ci 2 - 20a - 2c 2 . 



(49) 

- i„2 1 



In order to obtain the explicit form of the optimal control function from (49 1, by a variable replacement x = ±u 

/9-7C 



3 -, we rewrite this differential equation as (a.e.) 

(xf = 4x 3 - g 2 x - .93, 



(50) 



where g 2 = |[(/3 - 7Cl ) 2 + 3(/3 2 + 2 7 c 2 )] and g 3 = - 7C1) 2 - 9(/3 2 + 2 7 c 2 )]. 

From the classical theory of elliptic functions (see, e.g., [26, 27]), it is well known that the above differential equation 
is satisfied by the Weierstrass function g 2 ,gz) when the discriminant g\ — 27 'g 2 is nonzero. Therefore, one can 
express the candidate optimal control function as (a.e.) 



u{t) = ±4 



^(t + a;g 2 ,g 3 )- 



(51) 



1 (|M0) 2 



when f/2 — 27^3 ^ 0, where a = 

function u(t) turns at the point when u(t) cross the t axis. 

Consider the exceptional situations that the discriminant of (50) is zero, i.e., 

27 53 2 = (/3 2 + 2 7 c 2 )[(/3 2 + 2 7 c 2 ) -((3- 7Cl ) 2 ] 2 = 0, 
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92 , 53 ) ■ The sign of the above candidate optimal control 



(52) 



It is easy to see that either (i) (3 2 + 2jc 2 = or (ii) (j3 2 + 2^c 2 ) - (/3 - 7C1) 2 = 0. 

For the case (i), a further use of (44) leads to that 1 + 2jc 2 — (jua — f3) 2 — ^uc 2 = 0. Thus, if 7 < 0, we have 
Ua — ^, ub = ub(0) and uc = 0, which determines the candidate optimal control function by u(t) — ub(0) (a.e.). If 
7 > 0, from (49), we have (a.e.) 



7r 2 

4 [U 



+ -(/3- 7Cl )] 2 , 



(53) 



whose solutions can be expressed as follows. 
• If (3 — 7C1 < 0, then (a.e.) 



u(t) 



'2( 7Cl 


-P) 


7 


/ 2( 7 c 1 


-0) 



tanh 
coth 



V^P^ + a) , when |« fl (0)| < 



yf^Et + a), when \u B (0)\ > 



'2( 7 ci 


-P) 


7 


'2( 7 ci 


-P) 



(54) 



where 



Ar cth ( ± 



2( 7 ci-/3) 



«b(0) 



'2( 7 ci 


-0) 


7 


/2( 7Cl 


-0) 



(55) 



• If (3 — 7C1 = 0, then (a.e.) 



u(t) = ± 



Vl{t + a) • 



(56) 



where a = ±^ub(0). 
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If j3 — 7C1 > 0, then (a.e.) 



U(t) = ±4 



' 2(701 - d) 
1 



tan 



7C1 - /3 



t + a 



(57) 



where a = arctan 



(±v/^j^(o) 



For the case (ii), the equation (49 1 can be simplified as (a.e.) 

(uf = V + i(/3- 7 ci)h 
4 7 



(58) 



If 7 < 0, clearly, the above differential equation has no nontrivial solution other then u{t) = when j3 — 7C1 < 0. 
When (3 — 7C1 > 0, one can immediately obtain the optimal control function from (58 1 as (a.e.) 



«(*) 



ch (±V/3 - Tci* + a) ' 



(59) 



where a = Ar ch 
follows. 



ub(0) 



If 7 > 0, then from (58), one can obtain the candidate optimal control function as 



• If /3 — 7C1 < 0, then (a.e.) 

where a = arcsec ^| yT^zrpUB (0) 

• If /3 — 7C1 > 0, then (a.e.) 



tt(t) 



/7Ci ~ g , 
7 



(60) 



*(*) = 



sh (±V/3 - 7Ci< + a) 



(61) 



where a = Ar sh 



«b(0) 



Remark: In the case of a = and 7^0, the optimal control function can be expressed by the Weierstrass elliptic 
function only when the discriminant g\ — 27g 3 of (50) is nonzero. When g\ — 27g 3 = 0, the optimal control function 
is reduced to elementary functions. In comparison with the case that g\ — 27g\ 7^ 0, the case of g 2 — 27g\ = occurs 
with only a probability of zero. The similar arguments also can be made in the following two cases. 
2. Case a/0 and 7 = 0. 

In this case, control term B of system ([lj is parabolic. Accordingly, Eq.(48) can be simplified as (a.e.) 

{uf = au 3 + (3u 2 - 2acm - 2(3c x - 2c 2 , (62) 



With x = \au + jk(3, (621 can be recast as (a.e.) 



(i) 2 = 4x 3 - g 2 x - g 3 
where g 2 = ^(/3 2 + 6a 2 ci) and g 3 = ^(18a 2 Pd + 27a 2 c 2 - (3 3 ) 

9i 



(63) 



If the discriminant of (63 1 g\ — 27 'g\ ^ 0, one can again obtain the candidate optimal control function in terms of 



Weierstrass elliptic function (a.e.) 



u(t) 



1 



a 



4&(t + a; g 2l g 3 ) 







(64) 



where a = Sf" 1 (i[au B (0) + f ]; g 2 , g 3 ) ■ 

Otherwise, Ax 3 — g 2 x — g 3 has repeated zeros. Accordingly, the candidate optimal control function can be obtained 
from (62 1 or (63 1 as follows. 



9 



• If g 3 < 0, then 



u(t) 



'93 



l+cxp(± v /-6^g 3 "t + a) 
1 -cxp(±^/-6^t + a) 



where a = In 



• If (73 = 0, then 



where a = ±- 



V««b(o)+| 
• If ^3 > 0, then 



«(*) = 



4 g 
a(i±a) 2 3a' 



'S3 



3a ' 



3a ' 



(65) 



(66) 



(67) 



where a = arctan 



cm B (0)+§-4 3/g^ 



3. Case a ^ and 7^0. 

In this case, none of the drift term A and the control term B is parabolic. Let 

f(x) = ^2 4 + ax 3 + (P - 7Ci)x 2 - 2ac x x + 7C1 2 - 2/3ci - 2c 2 



■92 = ic+%ci + ^(/3-7 Cl ) 2 , 



53 



7(/3-7Ci) 
24 ' 



24 216 



22 9 
16 16 ' 



(68) 



where c = 7C1 2 — 2/3ci — 2c 2 . Let 27, 2 2 , 2 3 and 24 denote the roots of the equation f{x) = 0. 

When g\ — 27g| 7^ 0, it can be verified that Xi 7^ Xj (1 < i < j < 4). It is known that the solution of Eq.(49l still 



can be written down explicitly in terms of the Weierstrass function (see 26J chapter XX), which is given by (ale.) 

6f(xo 



u(t) = x + 



(69) 



24Sf(t + a; «? 2 , ff3 )- f"(x y 

where a = Sf" 1 ^ u f B ^°l Xo) + 37 f"{%a) 5 92, 53) and 2 € {21, 2 2 , 23, 24}. 

When — 27<7g = 0, the polynomial f(x) in (68 1 has repeated zeros. Then, there are four different situations 
accordingly. 

1) f(x) has only one 2-fold zero. 



In this case, without loss of generality, we can assume that 24 x\ = 2 2 7^ 23 7^ 24. Accordingly, (49) can be 
rewritten as (a.e.) 



(u) 2 = - (u — Xi) 2 [u 2 — (x 3 + x 4 )u + 2 3 2 4 ] 



From (70 1, one can compute the corresponding candidate optimal control function as (a.e.) 

27(2l - X 3 )(xi - Xi) 



u{t) 



|7(a; 3 - 2 4 )|sin (±5^-7(21 - x 3 )(x 1 - 24) t + aj - 7(221 - 2 3 - 24) 
where a = arcsin ll 2 ^^^ whcn l{xi _ x ^ Xl _ ^ < . and 

7(2xi-X3-a;4)+2( cxp(± A -v/7(:ri -rr 3 ) (xi -x 4 ) t+a) — 7< 2x i "^3 ~ x 4> ) -»/7(:e 1—2:3 )(:ri— 2:4) 

u(t) = ) V^i-sK-i-^ + Slj 

(cxp(±^ 7 ( 3 ; 1 - 3! 3)( a :i- 3! 4) t+a)- t("i-3-*4> _ 7 



(70) 



(71) 



(72) 
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where a = In 



y/'y(uB(0)-xi) 2 +-f(2x 1 -x 3 ~X4)(uB(0)-xi)+i(x 1 -x 3 )(xi-x i )+^/^(x 1 -x 3 )(x 1 -X4) 

ub (0)— xi 



y(2x 1 ~x 3 -x 4 ) 



2y -y(xi-x 3 )(x 1 -Xi) 

when "f(xi — x%){xx — x 4 ) > 0. 

2) f(x) has two different 2-fold zeros. 

In this case, without loss of generality, it can be assumed that x\ = X2 ^ X3 = X4. Accordingly, (49 1 can be 
rewritten as (a.e.) 



{uf = -{u - xx) 2 {u - x 3 ) 2 . 



If 7 > (there is no nontrivial solution for 7 < 0), from ( |73[ ), we have 

(xi — x 3 ) exp (±^(xi — x 3 )t + a) 



«(*) 



exp (±Z(xi — Xs)t + a) — 1 



where a = In UB ^~ Xl . 

u B (0)-x 3 

3) f(x) has one 3-fold zero. 



It can be assumed, accordingly, that x± = x 2 = x 3 x 4 . Then, (49) can be rewritten as (a.e.) 



{uf = -{u - x 1 ) 3 (u - x A ). 
The corresponding candidate optimal control function is given by (a.e.) 

47(0:1 - Xi) 



u(t) 



1 2 {x 1 -x i f{a±Uf-4-f 



Xi, 



where a = - l{xi _ Xi) 2 (ub(Q) _ Xl) VtMO) - xi) 2 + Wi ~ Xi){u b {0) - si). 
4) f( x ) has one 4-fold zero. 



In this case, we can assume that x\ = x% = X3 = X4. Then, (49 1 can be rewritten as (a.e.) 



{uf = l{u- Xl )\ 



(73) 



(74) 



(75) 



(76) 



(77) 



If 7 > (there is no real solution for 7 < 0), from (77 1 one can obtain the candidate optimal control function as (a.e.) 

1 



u(t) 



a±^t 



+ Xi, 



(78) 



where a 



i(0)-xi 



V. EXAMPLES 

In this section, we give two examples for illustration. 

Example 1: Consider the case when A = K x + 2K Zl B — K x and the target evolution matrix has the form 
X f = e BK * (6 e R). 

It can be checked that the abnormal optimal control it a b norma i(£) = — |g fit | = — 1 steers system (1 1 from the initial 
state I2 to the final state Xf in time Tt — (6 mod 4tt)/2, with the performance measure given by 

1.(0 mod 4tt)/2 



J (^abnormal) 



Wabnonnal(*)^ = (0 mod 4tt)/2. 



(79) 



Actually, taking the matrix S as 



S = s z K z , 



(80) 
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we have 



H(S; A o5 u; X {t)) = (S, X^i)" 1 [A + u{t)B]X\t)) 
= (s z K z ,e- 2tK » [K x + 2K Z + uK x ]e 2tK *) 
= (s z K z , (1 + u)[cos(2t)K x - sm(2t)K y ] + 2K Z 
= 2s z . 



(81) 



which is minimized with respect to u by the abnormal optimal control M a bnormai(i)- 

Example 2: Suppose that the system Q is given by A — K z and B = — K x + K y . Consider the optimal steering 
problem with the terminal condition Xf = e~ 2Kx+2K v _ 
Write the matrix S as 



S S X K X + SyKy + s z K s 



(82) 



In order to achieve target evolution matrix optimally, according to the Eqs. (46 1 and (47 1, the coefficients s x , s y and 
s z in (82 1 are required to satisfy 



sinh(2^) 
Sz ~ v/2[cosh(2\/2) - 1] {Sx + Sy '' 



(83) 



Fig [T] shows the resulting evolutions. 
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10 20 30 , 40 50 60 



FIG. 1: The Distance Between the State X(t) and the Target Xf at time t. A numerical treatment shows that when the 
two independent parameters (s x , s„) take the values of (-2.09895,0.99801), (-0.844738,0.406815), (-0.523766,0.254595), 
(-0.379568,0.185591), (-0.297326,0.145922), • • • , the target evolution matrix X f = e ~ 2K *+ 2K y can be reached at time Ti = 

9.625, T 2 = 21.950, T 3 = 34.395, T 4 = 46.880, T 5 = 59.394, • • ■ . Where the norm \\X\\ F is defined as \\X\\ F = \ x a\ 2 - 



Correspondingly, the control functions given by 



|ui(t)| = JW{t + 0.5740 - 1.2502z; 9.8362, 4.4871) + 2.4990 

_ 1.0547 (84) 

• V /S?(t+0.5740;9.8362,4.4871) + 1.2495' 



\u 2 (t)\ = J2&(t + 0.6705 - 2.0412i; 3.6531, 1.3108) + 1.2418 

_ 0.3487 (85) 

• v /»(t+0.6705;3.6531,1.3108)+0.6209' 

\u 3 (t)\ = ^2<g(t + 0.7051 - 2.522H; 2.6390, 0.8203) + 0.9956 

_ 0.2045 (86) 

^/S?(t+0. 7051;2. 6390, 0.8203)+0. 4978' 
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|u 4 (*)| = y/2&(t + 0.7247 - 2.8643i; 2.2391, 0.6435) + 0.8954 



y/ {?(t+0. 7247:2. 2391, 0.6435)+0. 4477' 



(87) 



|u 5 (i)| = yf2&{t + 0.7370 - 3.1118z; 2.0254, 0.5543) + 0.8416 

_ 0.1115 

^W(t+0. 7370;2. 0254, 0.5543)+0. 4208' 



as shown in Fig [2] are candidate optima. 




10 20 30 ( 40 50 60 70 
FIG. 2: The Candidate Optimal Control Functions. 



The corresponding performance measures are shown in Fig [3| It is easy to observe that there is a tradeoff between 
the consumed time T and the cost index J (it). 
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FIG. 3: The Performance Measures Corresponding to the Candidate Optimal Controls. Where J(ui) = 7.3473, J{u?) = 3.1318, 
J {us) = 2.0125, J(ui) = 1.4932, J {us) = 1.1873, ■ ■ ■ . 
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According to the decomposition algorithm given in [TJ], piecewise constant control laws can be designed to achieve 
the target evolution Xf = e ~ 2K ^+ 2K v gg well. One possible control law is given by 



v 



0, when t E [0,ii + 2ni7r); 
( t )= { 75' whent G [ti+2n 1 7r,ti+t 2 + 2ni7r]; (89) 
0, when t S (t x + t 2 + 2ni7T, 2t x + t 2 + 2(m + n 2 )ir], 



where ni,n 2 £ N + , c > 1 and 



<i = 2arccot (-ccth V2 - \J c 2 cth 2 V2 - 1 
f 2 = -r=r Ar cth f -4=r Vc 2 cth 2 V2 - 1 



The corresponding performance measure is given by 



■2 



J(u) = —== Ar cth - Vc 2 cth 2 \/2-l . (90) 

It can be verified that J(u) increases monotonously with the increase of c. Since 

lim J(u) = 1 » 1.9351, (91) 

c ^ Vcth 2 V2 - 1 

as a comparison, the performance measures J{u±) and J(u$) are approximately 13 and 39 percent, respectively, less 
than J(u). 



VI. CONCLUSION 



In order to minimize the decoherence effect, energy optimal control problem for the quantum systems evolving on 
the noncompact Lie group S^l, 1) is taken into account in this paper. We showed that explicit expressions for the 
optimal control functions can be obtained analytically. To minimize the considered quadratic performance measure, 
the control functions with respect to abnormal extremals are constant functions of time t, while those with respect 
to normal extremals are expressed by the Weierstrass elliptic function. 
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